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abstract 


An on-line spatial light modulator developed at CMU suitable as the input 
transducer for a real-time optical data processing system *- s .described 
(Chapter 1). 1000 x 1000 point resolution and on-line operation at TV frame 

rates is possible. The use of such a device in the analysis and processing of 
radar data in real time has been considered and initial results reported. 

The fundamentals of signal processing for linear (Chapter 3) and planar 
(Chapter 4) phased array radar data as well as Pulsed Doppler and Linear-Step 
FM system (Chapter 7) using this -input modulator have been developed. Time 
and space multiplexing input formats, full and thinned aperture processing, 
and amplitude and phase modulation schemes using the system are developed, 
analyzed and compared. 

An interface from the optical processor to a control-digital computer 
has been designed, constructed, and tested. The ability to digitize the entire 
Fourier transform plane or the correlation plane of an optical processor in 
real time and input the resultant data to a digital computer for storage or 
for use in controlling the entire optical processor has been demonstrated. 

The cell resolution and size and the gray scale levels in the backplane are 
variable urtder program control. 

Real-time processing of linear phased array under data has been performed 
for the first time by other than an acousto-optic modulator. The Pockels 
tube modulator offers advantages in fabrication, angular coverage and resolu- 
tion, in the larger array size it is now possible to handle and in-depth of 
modulation etc. 

The input transducer, optical system, and computer interface have all been 
operated in real-time with real-time radar data with the input data returns 
recorded on the input crystal, processed by the optical system and the output 


plane pattern digitized, thresholded and outputted to a display and storage’ 
in the computer memory all in real-time synchronized to TV frame rates. 

No other input modulator offers the real-time operating characteristics 
of the Pockels tube. Its ability to process radar data in real-time has been 
demonstrated.' The correlation with theory and experimental results have been 
most encouraging. Recommended directions for future research in the area of 
optical processing of radar data in real-time are sketched in Chapter 7 and the 
proposal. 


CHAPTER 1 


OPTICAL PROCESSOR 
1.1 Introduction 

This contract involved a one man year effort to investigate the feasi- 
bility of applying optical data processing techniques to the analysis of 
radar data. A theoretical analysis of the required data processing is inclu- 
ded; input formats compatible with an electron beam addressed input modulator 
are analyzed and presented; an interface from the optical processor to a 
general purpose minicomputer was constructed and is described; a phased array 
simulator was designed and built to serve as an input data source and finally 
experimental results of real-time radar data processed electro-optically and 
interfaced to a control minicomputer were obtained and are also included 
along with recommendations for future work in the area of real-time processing 
of radar data using optical techniques. Radar Data Processing in real-time 
was shown to be possible using C-MU's (Carnegie-Mellon University's) hybrid 
Pockels tube real-time spatial light modulator (See Section 1.2). 

As outlined in the proposal, the scope of the project included six areas: 

A. Formal extension of the theory of two-dimensional phased array radar pro- 
cessing with an electron-beam addressed light modulator. 

B. Determination of the ideal system configuration by means of evaluating 
simulations of various data formats. 

C. Automation of the necessary preprocessing of the electrical signals from 
a phased array antenna, and the post processing of the resultant diffrac- 
tion pattern. 

D. Design of the input logic network which interfaces to the video input of 
the electron gun. 

E. Investigation of feedback control loops to control the input plane format 
by sampling, analysis, and processing of output plane data. 
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Construction of peripheral electronics which convert investigators opt 
cal data processor to the task of a radar data processor. 
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1.2 ON-LINE INPUT TRANSDUCER 

Optical data processing in real-time is nearing reality. One of the 
major obstacles in the implementation of such a system has been a two- 
dimensional electrical-to-optical input transducer from the real world 
to the optical processor. This device must be capable of creating a 
two dimensional image of an input electrical signal. A raster scan format 
and an input video signal are the easiest situation to envision. In this 
case a transparent version of the transmitted video image is formed. The 
material on which this image is formed must be capable of spatially modu- 
lating the amplitude and/or phase a collimated input laser beam. 

Figure 1.1.1 shows a block diagram of the input transducer used. The 
input electrical signal, assumed for simplicity to be of TV raster scan 
video format, modulates the cathode and hence the beam current of an off- 
axis electron gun. The deflection of the electron beam is synchronized 
to horizontal and vertical sync pulses. The target is not a phosphor, 
but rather an electro-optic target crystal, on which the modulated electron 
beam deposits a variable charge pattern. A transparent conducting layer 
on the back surface of the crystal serves as the system's ground reference. 

The mechanism by which this target material modulates light has been 
well documented. The locally varying voltage across the crystal at any 
point is a controllable function of position. The electric field produced 
at every point causes a retardation between components of the light beam 
incident at that point. By the linear longitudinal electro-optic orPockels 
effect, this controlled voltage can be used to modulate the phase of the 
incident coherent light beam at each point on the crystal's surface. 

When placed between crossed polarizers this phase modulation can be 
converted to amplitude modulation. The input and exit windows allow the 
light beam to pass through the target in a transmission mode. 

•*r 

*r 



Fig. 1*1.1 Block Diagram of. Input Transducer (Pockels Tube) 
Code: 1 = optical window, 2 = write electron gun, 

3 = target crystal (KD 2 P) , 

4 = Conducting transparent grounded electrode. 



Electron beam technology limits the obtainable spot size at the re- 
quired beam currents to about 1 mil. Most data and signal processing appli- 
cations require 1000 line resolution. 1000 x 1000 point resolution is 
possible by Curie point operation of the target crystal. This high 
resolution is purchased at the price of storage of the charge on the 
crystal's surface. Real time operation necessitates the use of an erase 
gun which emits a defocussed spray of electrons toward the target at a 
low accelerating potential such that the crystal's secondary emission ratio 
is greater than 1. 

\ 

The detailed variation of the physical parameters of the target 
crystal with temperature and additional advantages of Curie Point operation 
of a Pockels tube are contained in Reference (2). Reference (3) explains 
the system more fully. The tube has been rather extensively documented 
elsewhere, References (l)-(6) , and thus only a brief explanation is included 
here for completeness. 

1.3 ELECTRO-OPTIC EFFECT 

The electro-optic effect is well known (7). Figure 1.3.1 shows the 

geometry of the electro-optic system. A 0° Z-cut thin basal plane of 

KDP of thickness L is shown with a voltage V applied in the direction of 

the optic axis z. The x and y axes are the crystallographic axes while the 

x' and y' axes are the electrically induced principle dielectric axes in 

the crystal due to an electric field E applied in the z direction. 

z 

1.3.1 AMPLITUDE MODULATION 

In crystals exhibiting the linear electro-optic effect, the indices 
of refraction for light polarized in the x' and y* directions are changed 
by an amount An in the presence of a longitudinal electric field 
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Output polarizer 
i ll to y ) 


Fig. 1.3.1, Amplitude Modulation Using the Pockels Effect. 


n ' = n - An 
x o 

n ' = n + An 

y o 


where n is the ordinary index refraction and n ' and n 1 are the indices 
o J x y 

of refraction along the x' and y' axes in the presence of an electric 

field. KDP and its isomorphs are biaxial in the absence of an applied 

voltage with n x = n^ = n^. The change in index of refraction is linearly 

related to the electric field E = V/L by 

z . J 

An = n ^r,,E /2 
o 63 z 


where r^ Is a tensor component of the index ellipsoid. 

This induced birefringence results in a differential phase delay or 
retardation between the components of an incident laser beam along the 
two induced axes given by 


3 

2Tm r c ,E L „ 
_ o 63 z V 

r = s * y- 


TT 


The phase retardation increases linearly with the voltage applied across the 



crystal until at V = V the phase retardation between component waves 

7T 

is ir and destructive interference can occur. 

The transmission through the polarizer Crystal/ analyzer system of 
Fig. 1.3.1 vs. the voltage applied across the crystal can now be evaluated. 

By combining the projections along the crystallographic axes when the 
polarization of the input beam is oriented along one of the crystallograph- 
ic axes at 45° to the induced electro optic axes. The amplitude transmission 
is then 

v 

A = sinS 

where S = T/2. The intensity transmitted at each point is 

I = I Q sin 2 6 

where 1^ is the incident intensity and 6 is the modulation which is 
linearly proportional to the applied voltage. 

1.3.2. PHASE MODULATION 

The electro-optic crystal can also be used to phase modulate the 

incident light beam point-by-point. Figure 1.3 .2 shows the polarization 

\ 

\ 

setup required for this scheme. The incident light is now linearly 
polarized along one of the induced electro-optic axes rather than along 
the crystallographic axes as in the amplitude modulation scheme of Fig. 

1 . 2 . 

The phase retardation along the x' axis is related to the electric 
field E z = V/L by 



Fig. 1.3 .2 Phase Modulation Using the Pockels Effect. 


where the terms are as previously defined. In Chapters 3 and 6 these two 
modulation schemes, amplitude modulation and phase modulation will be 
evaluated as they apply to optical processing of radar data. 


1.4 Extent of Research Performed 

i 

'■ " The lack of any available radar data, phased array or other, from 

i 

! NASA severely limited the scope of the work and the impact of the results 

i 

1 obtained. Any preprocessing needed would depend greatly on the specific 

radar used. Tradeoffs exist between the data availble to us, the format 
in which data could be obtained given complete access to a specific radar 
system, and the time and cost involved in constructing a hardware pre- 
processor. This preprocessor might be both special purpose (often for 
one radar system that might not be of adequate interest) and would in 

essence amount to converting data from one form (e.g. paper tape) back 
i 

into its original form (e.g. an electrical signal vs. time). These 
factors and the interest of NASA in phased array radar processing specifical- 
ly clearly indicated that no preprocessor be constructed. Rather a linear 


■ t . 
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phased-array waveform generator/simulator (Section 6.2) was designed, 
constructed, and used as the input data source. 

Only in this way could actual real time radar data processing be 
demonstrated. No contract renewal was indicated and hence all work was 
organized to present a unified treatment of the theory of 1-D and 2-D 
phased array radar processing, an evaluation of various input data formats 
for an electron beam addressed spatial light modulator, an extension show- 
ing the possibilities of such processing, and an actual demonstration of 
real time on-line radar processing using optical processing. 

i 

Chapter 2 provides the theory of beam steering, the optical genera- 
tion of such systems, and the analogy between beam steering and beam 
forming. Chapter 3 extends this treatment to the processing of linear 
phased array antenna data. Various input formats, modulation schemes, and 
the expected output data formats are theoretically treated. In Chapter 
4, the theory for a 2-D phased array is developed. Treatments of various 
input formats, optical configurations, and output patterns are again 
included. Various new theories such as a scheme for simultaneously ob- 
taining azimuth and doppler data are included. The digital computer to 
optical computer interface is described in Chapter 5. In Chapter 6, 
various phase and amplitude modulation schemes are included. On-line real 
time linear phased-array radar processing results are shown, performed 
optically for the first time by other than acousto-optic techniques. The 
interaction in real time of the output optical/digital interface is also 
shown by showing images of the on-line digitized versions of the output 
patterns as stored and digested by the computer. 

Possible extensions of this work were submitted separately to 
Reginald Inman and Joe Kerr of MSFC as an informal proposal for continued 

■tt- 


rt . ■ 
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real-time optical processing of radar data. Topics proposed include: 

1. Pulsed doppler radar processing in real time by optical techniques 

2. Linear FM stepped radar processing in real time by optical techniques 

3. Analysis of computer interface requirements 

4. Spatial filter study for real time correlation in radar processing 

5. Multiplexing 1-D transforms for radar processing 

6. Automatic real time range gating by optical correlation 

7. Real time fine range resolution by optical correlation 

8. Improved S/N ratio schemes using optical correlation 

t 

9. Real time CW dopper processing using optical techniques 

10. Simultaneous azimuth and doppler data from one processing frame 

11. Alternate electron optic scanning scheme, optimized for radar processing 

We now have access to chirp and doppler data and would hope to see 
these proposed efforts, all. of which are proposed and currently possible 
in real-time, pursued since we are standing on the edge of major demon- 
strations in this area. 



CHAPTER 2 


BEAM STEERING IN PHASED ARRAY ANTENNAS AND OPTICAL BEAM STEERING 
2.1 Array Antenna Characteristics 

Essentially, array antennas consist of a number of indi- 

vidual radiating antenna elements, spaced in accordance with some given 
pattern. The characteristics of such antennas are determined by the 
geommetric position of the radiators. The amplitude and phase of the 
excitation are modulated by controlling the individual elements to give 
beams that are scanned electronically. Three different methods for 
scanning are possible: phase shift, frequency shift and time delay. 

We shall concentrate our attention on phase shift scanning and show the 
analogy between such beam steering and optical beam steering using 
"optical" antennas. This analogy is the basis for the optical proces- 
sing of signals received at phased array antennas, or beam forming, 
via the reciprocity theorem. 

Linear as well as two-dimensional antennas will be considered. 

For a uniform linear phased array of N radiating isotropic elements, 
let the transmitted signal at the nth array element be of the form 


where 


v n (t,<}>) = p T (t)cos 2 7T [ f c t .- n<j, o ] 


1, |t| < T/2 

- , . • where T is the transmitted pulse width 

0, otherwise r 

f c = RF carrier frequency 

♦ = differential phase delay in RF cycles applied per element 


The entire signal transmitted by the array is thus a waveform whose equi- 
phase front forms an angle 0^ with the normal to the array — or boresight 


direction — where 
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4> X 

sln0 = — — 
o d 

X = free soace wavelength of KF carrier 
c 

d = equidistance spacing between array elements 
Fig. 2.1.1 describes the scanned array configuration. For convenience 
in the computation, the signals v (t) will be taken in its analytic form 


V n (t ^o ) = p T (t)e 


j 2 tt [f c t - nct> o ] 


with 


- Re{V n (t,«| )n )} 


Application of the Fresnel-Kirchoff scalar diffraction theory to the antenna 

shows that, in its far field (this is for space points at a distance great- 
2 2 

er than 2N d A c ) » the signal received by a point located on a direction 

forming an angle 0 with the boresight axis is 

N-l j2r[f (t - t ) - n<f> + n<j> ] 

V(t P T (t - t )e C K 

n=0 


j2irf (t - t„) N-l - j 2 un (<p - <f>) 

= p T (t - t R )e C R I e 


n=0 


where <#> — (d/X c )sin0 , and t is the one-way target range delay 

- R 

By adjusting the phase terms relative to the center of the aperture, 
we obtain 


V(t, *,* ) = <Kt)E (0,0 ) 

O 3 0 


where ip(t) = p T (t - t R )e 
and 


j 2^rf ( t - O 


N-l 

2 -j2im(<{> - <{>) -sinNTi((|> - tp ) 


2 


sin[Nir(d/A ) (sin0 - sin.0 )] 

= 2__ 

sin[n(d/A ) (sin0 - sin0 )] 
c o 

For 6 q = 0, the expression 

sin[N7r(d/X )sin0] 

T? / Q \ = £ 

a v sin[ir(d/X c )sin0} 


is known as the array factor for isotropic radiators. When the radiating 
elements are not isotropic but have a radiation pattern E e (0) , the antenna 
radiation pattern is given by 


E(0) = E (0)E (0) = E (0) 
e a e 


sin[NTr(d/X )sin0] 
c 

sin[ir(d/X )sin0] 
c 


Examination of the array factor shows the following points* 

1) E fl (0) * s repetitive, with grating lobes spaced at values of 0 

for which Tr(d/X c )sin0 = nir, n = 1, 2, 3, ... 

*1 

2) The half-power (3 dB) lobe width has the value ^ thus application 

of. the Rayleigh resolution criterion yields N resolvable main lobe 
positions corresponding to N values of sin0 for 0 between 0 and tt radians 

Fig. 2.1.2 illustrates the former results. An approximation to the 
array factor for small values of 0 is given by 


sin[TrN(d/X) sin0 ] 

E(e) - 

uN(d/X)sin0 

c 


which corresponds to the Fourier transform of the illumination across the 
continuous aperture. This approximation holds when the array aperture is 
larger than several wavelengths. For larger values of 0 the obliquity 
factor must be included so that 
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. . . 1 + cos0 . sin[uN(d/X)sin0] 

E < 0) ~ < 2 } uN (d/X^ sinQ 


For the scanned array system, the array pattern was shown to be 


sin[Nn(d/X c ) (sin0 - sin0 Q ) ] 
E a^ 0 ^ ” sin[ir(d/X ) (sin0 - sin0 )] 


whose grating lobe maxima repeat at values of 0 such that 


ir v( s i n 6 “ sin0 ) = mr, n = 1, 2, 3,..,. 

X c ° 

The 3 dB antenna beam width is again ^ y ^ . A total of N values of sin0 Q 
can thus be resolved between two consecutive grating lobes of the array 
pattern. Fig. 2.1.3 shows this result. 

Again, the antenna pattern can be approximated by the Fourier transform of 
the illumination across the continuous aperture multiplied by the obliquity 
factor 



E (0) = ( 


1 + COS0 
2 


) 


sin[irN(d/X ) (sin0 - sin0 )] 
c o 

nN(d/X ) (sin0 - sin0 ) 
c o 


and the power associated with this pattern is 


|e (9)i 2 - ( 1 * ;° s9 > 


2 sin [TTN(d/X ) (sine - sin0 )] 
c o 

TT^N^(d/X )^(sin0 - sin0 ) ^ 
c o 


For 0 <_ 30°, however, the term ( — + ^ c - ° - s ~- ) has a maximum value of 1 
at 0=0 and decreases smoothly to 0.93 for 0 = 30° at which 

o 2 

( 2 ) = 0.87. The power pattern for near-to-boresight direc- 

tions is thus approximately 

2 sin^[N7r(d/X ) (sin0 - sin0 )] 

|E a (e)| 5 f— 

[irN(d/X c ) (sin0 - sin0 Q ) ] 

where we recall. that 0 q corresponds to the steered direction of the antenna 
and 0 to the angle with the boresight axis as shown in Fig. 2.1.1 
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From this analysis it is clear that the ideal number of resolvable 
beam positions in the far field on the antenna is related to the array 
factor which in turn is determined by the antenna geometry. This array 
factor is thus a measure of the maximum angular resolution from an ideal 
processor. Once this spacing d between elements and the carrier wavelength 
Xc are chosen the maximum sector scanning angle 0 for non ambiguous 
positions is defined by (d/X^) sin 0 <_ ir. Within + 0 the number of 

resolvable beam positions is N, the length of the array; e.g. for d/A^ 

< 0.53 N beam positions can be resolved between 0=£6O°, while for d/A c ' 

<0.5 yields N resolvable positions for sin 0 with 0=+ 90°. It will be assumed 
from now on that the relation d/X £ = 0.5 holds for all arrays considered. 
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2 . 2 Signal Processing-Requirements 

Phased array radar antennas present the advantage^ over conventional 
antennas, of rapid beam steering, this being achieved by varying the 
relative phases of the excitation applied to the different elements in 
a precise manner; no inertia, friction or moving parts are source of 
problems with this mode of operation. 

With the capability of rapidly and accurately positioning phased 
array antennas, it has become possible to develop radars that perform 
multiple functions interlaced in time. Sequential scanning in azimuth 
and elevation can produce the scene in front of the antenna as well as 
range and velocity information. 

Another mode of operation is also possible with a phased array that 
can only be realized with this type of antenna. The whole solid angle of 
space to be explored can be illuminated at once by using a broadbeam 
antenna as transmitter, and' the desired information about targets in the 
searched volume can be extracted by processing the signals received at 
every element of the array. Elevation angle can be obtained for linear 
arrays, and elevation and azimuth angles for planar arrays. The require- 
ments placed on the processor are, however, quite formidable: real time 

operation for a broad class of targets require parallel processing of all 
array elements as well as multitarget capability. 

Essentially, the function to be performed by the processor, for 
every range bin of interest, (a linear array is considered), consists of 
combining the signals received at every element of the array to obtain an 
output signal V(.9) such that for any angle 0^ at which a target is present 
the output presents a maximum and for all other values 0 4 0^ it varies 
in proportion to the antenna radiation pattern associated with the mode of 
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operation of the antenna, (pulsed RF carrier in most cases) . The processor 
thus must correlate the phase slope present at the returned signals with 
all possible linear phase slopes and determine, for every peak of the correla- 
tion function, the value of the elevation angle of the target. This can 
be accomplished by 

1) Quantizing the angle coverage in N discrete values 0^* • • 8^ 

2) For every such value 9^, delaying each of the RF signals by an 

amount such that the nth signal is delayed by an amount nt , 

where T = (d/X ) sin 0. 

c 1 

3) Sum the resulting ensemble of N signals which are now of the 
same phase 

4) Coherently integrate the summed signal for the entire signal 
duration. 

It can be seen that the processor must perform the aforementioned 
operations for all values of t of interest, i.e. for all values of the 
angle variable 0. 

In Section 2.1 the resolution capability of a phased array was derived and 
its relationship to the array geometry demonstrated using the array factor. This 
process is beam steering. We now consider the analysis required on the signals 
received at the detector array from a- point target. This processing is beam 
forming since the required analysis involves shifting, delaying, and summing 
all returns to form all possible beam directions in which targets may exist. 

There clearly exists a duality between beam forming at the detector array and 
beam steering at the transmitting antenna array. This reciprocity will become 
apparent from a detailed discussion of an ideal phased array processor. 

To describe the process more fully let the- line antenna of length D 
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consisting of N elements spaced at distances d, with illumination function 

I(x) e j2lr x (x) 

where I(x) accounts for the possibility of weighting the signals received 
at every antenna element and x(x) is the phase delay (in RF cycles) intro- 
duced at the receiver. The variable x quantized at discrete values 

corresponding to the locations of the N elements 

N-l 

x - id, i = 0, + . . . + —• 

Let the antenna receive a plane wave impinging under the angle 0, as 

v 

depicted in Fig. 2.2.1. The processor output signal is the sum of all 
contributions, weighted by the illumination function and delayed by amounts 
depending on the path length differences. The differential delay for each 


,j*- v 


r * ■ 
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element is 


t = — sin 0 
c 


In the general case of a narrow band signal with complex envelope 
y(t) , the signal received at a particular point on the antenna has the 
form y (t-x)e^ 7T ^c^ t T \ The processor output signal is then 


V(t ,0) = l 


n =_N-1 
2 


y(t- — sin 0) e 
c 


j 2 tt ( f c t sin 0)i(nd)e' i27r X ' Rd) 


For waveforms with small percentage bandwidth, the complex envelope , v 

distortion is negligible across the array elements. This allows us to write 

N-l 

V(t,0, 0 ) - p(t) E (0) 2 l e [sln ». - sln « 

n = N-l 
"2 

where the element pattern has been included to take into account the 
directivity of each receiving element, and where the following assumptions 
have been made: 

1) I(nd) = 1 f or all n: that is, uniform amplitude illumination across 
the entire array . 

2) )^nd) = n~ sin 0 q : that is, the incoming signal is correlated 

c 

against a linear phase slope corresponding to an assumed incidence angle 0 q . 
The summation term in V(t, 0, 0 q ) reduces to 

sin [N ttOt— ) (sin 0 - sin 0 )] 

o 

sin[Tr(-T--) (sin .0 - sin 0 ) ] 
c 0 

which corresponds exactly to the scanned array pattern derived in section 
2.1. This part of the output of the ideal processor has thus been shown 


to correspond to the scanned array pattern of the receiving antenna. The 
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Fcurier transform relationship between the approximate array pattern and 
the aperture illumination function of the antenna shall be conveniently 
exploited in the implementation of the optical processor. The same 
conclusion can be arrived at by writing V(t ,(f> ,<j> o ) as 


N-1 J 

VCt.M ) - P(t)e 32,rf c t (^£SS4) / 2 I( X )e 32 " x(x) e- 32 * x sin< * dx 

° N-l , 


where x = r— ■ = (— ) - , i = distance across the array, the summation over a 
a c t 

c c 

discrete set of points has been approximated by an integral and the element 

14-cos 0 

pattern equal to the obliquity factor — z . The antenna pattern is then 


E (0, (f ) - (it£2=2) / " I(x) e 32ir * (x) e - j2nx sln0 dx 

o l * 

—oo 

1 | | <D £sin0 

which, for I(x) = L* ~2. and X (x) = — r = x sin 0 (where D = 

* u, otherwise X o 

c 

(N-l)d) reduces to: 


E (0,0 O > - f P D (X) e 3 2 *>‘(s 1 ''» 0 - » dx 

— CO 

and thus, 

V(t,0, 0 ) = y(t) e^ 2 "^ E (0, 0 ) 
o o 

where a Fourier transform relationship holds between E (0, 0 Q ) , the 
approximate scanned array pattern, and the illumination function of the 


antenna . 
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The approximation of the element pattern E e (equal to /cos <j> for a well- 
designed array) by the obiquity factor — is valid for values of if> up 

to 60° or 70°. 

. The Fourier transform relationship which was just proven to exist between 

i 2 tty (x) 

the illuminating function I(x)e J A and the antenna pattern E (<|>,<{> o ) is the 
basis for optical beam steering which we will now consider. The analogy be- 
tween beam forming and optical beam steering will then be apparent. 



/ 
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2 . 3 Optical Beam Steering 

It is well known that optically a lens forms in its back focal plane 
a light distribution — in amplitude and phase-which equals the Fourier trans- 
form of the light distribution existing in its front focal plane, when mono- 
chromatic coherent light is used. In the physical situation, if a collimated 
plane light wave is modulated with the one dimensional phase distribution 
corresponding to the linear phase slope of a RF waveform impinging over a 
linear array antenna and then an one-dimensional Fourier transform is taken, 
the output plane presents a light amplitude distribution which corresponds to 
the continuous aperture antenna pattern (for appropriate scaling factors) and 
which presents a maximum at a location corresponding to the angle of incidence 
on the array of such waveform. Fig. 2.3.1 depicts this situation. 

Also, for a planar phased array, if the two-dimensional phase distri- 
bution of a plane waveform impinging on the array at angles 0 and as 
described in Fig. 2.3.2 is used to modulate the phase of a plane collimated 
light wave by using a convenient transducer, the Fourier transform of such 
light distribution is again the continuous aperture antenna pattern which 
presents a maximum at a location corresponding to such incidence angles 0 and 
ip as shown in Fig. 2.3.3. 

Optical beam steering based on such principles has been successfully 
performed with devices such as the Membrane Light Modulator at A.P.L. 

(Johns Hopkins University). (The reader is referred to Report TG 1193A for 
further details). This mode of operation requires, however, knowledge of the 
phase of the RF signal received at every array element prior to the operation 
of beam forming. For an array with N elements, a total of N phase detectors 
and D/A converters are required for parallel operation, and there is no 
provision for optical waveform processing, since the input transducer is a pure phase 


Linear phase 

slope input Cylindrical lens Spherical lens 



Fig. 2.3.1 Optical Beam Steering for a Linear Phased Array Antenna 
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modulator. It seems thus advantageous to build a processor which operates on 
the received waveforms themselves rather than on the RF phase. Most of the 
pre-processing hardware is eliminated and the required operations of amplitude 
weighting and filtering are much simpler. ' 

With the availability of the Pockels light modulator, the possibility of 
such processing exists. This has been the center of much of our recent research 

I 

effort. 
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CHAPTER 3 

ELECTROOPTICAL PROCESSING OF LINEAR PHASED ARRAY ANTENNA SIGNALS 
3.0. Introduction 

In Chapter 2, the Fourier transform relation between the continuous aper- 
ture antenna illuminating function and the far-field antenna pattern in beam 
steering was shown to be the reciprocal of the relation between the array illum- 
ination function of the receiving antenna and the processor output in beam form- 
ing. One possible processing format which related beam forming and optical 
beam steering was also presented. In this chapter, beam forming by an optical 
processor is shown to result when the returned electrical signals, converted 
to optical signals by means of an appropriate transducer, are used in optical 
beam steering. It will be shown that such optical beam steering yields an out- 
put corresponding to that of the ideal processor defined in Chapter 2. The op- 
tical system, several possible input formats to be used with the input trans- 
ducer described in Section 1.2, and derivations of the format the optical output 
pattern are also included. 
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3. 1 Signal Preprocessing 

The simplest configuration for a linear receiving array was shown in Fig. 
2.2.1. The RF signal reflected from an assumed point target located at a 
direction forming an angle <j> with the boresight axis at a range R is repre- 
sented by a plane waveform impinging on the linear array. The elements of the 
array are equi-spaced by distance sd. The ideal processor described in Chapter 
2 is shown on the same diagram; it consists of N delay lines, a summer and a 
coherent integrator yielding the processor output V ( <j>) . 

Modern radars transmit RF signals of frequencies in the S, C and X bands. 

\ 

With the limited input plane resolution available and the present electron gun 
scan at video rates, it is impossible to directly record the returned RF sig- 
nals, The approach usually taken is to first amplify the RF returns in a radio- 
frequency amplifier, then heterodyne down such signals with a local oscillator; 
the resultant IF signals are range gated before processing to extract angular 
information. Fig. 3.1.1 shows a block diagram of such a system. 

3.2 Input Format 

The most feasible input format is now described, with alternate configura- 
tions to be discussed later. To process linear phased array antenna data, the 
time history of each detector/receiver element after amplifying and heterodying 
is recorded along each horizontal line in the input plane. With a total of N 
detector elements, N horizontal lines are used on the input plane electrooptic 
crystal. The spacing between these lines is proportional 
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to the spacing of the detector elements in the detector array. Such an input 
format is depicted in Fig. 3.2.1. 

Let us now consider this processing and the resulting output plane 
pattern in detail. Recall that the input electrical signals from the 
detector elements will be used to modulate the beam current in an off- 
axis electron gun. This gun then deposits a charge pattern on the electro- 
optic target crystal which results in a voltage being developed across 
the crystal and hence an electric field all of which are controllable 
functions, point-by-point, of the input electrical signal. By the linear 
longitudinal electro-optic effect, this charge distribution then modulates 
the amplitude amd/or phase of a collimated input laser beam, as described 
in Section 1.3. 

For an array consisting of N identical elements, the optical image placed 
on line n of the crystal should be a representation of the (heterodyned) 
signal received by the nth antenna element v in C fc ) over the integration time 
T. for. the range delay of interest. The integration time T. is de- 
fined as 



where 

= time at which recording of returned signals begins 
= time at which recording of returned signals ends. 

With the electron gun scanning the electrooptic crystal at horizontal TV 

rates, T. is approximately 50 psec. Fig. 3.2.2 shows the pulse width T . 

1 P 

range delay t and integration time T. for one detector element and their in- 
R i 

ter-relations. The value of T^ is set by the range gates following the IF 


mixers in Fig. 3.1.1. 




t=o t=0 


Fig. 3.2.2 Relationships Between Integration Time.T. , Range Delay t 

1 R 

and Transmitted Pulse Width T . 

P 
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Let the signal transmitted by the antenna be an RF pulse of the form 

v T (t) = p T (t) cos (2irf ^t) 


where 


P T (t) = 


l,|t| < T p /2 
0, 1 1 j > T p /2 


After being reflected by a target at a range R, the received signal at the nth 
element of the array is given by 


= p T (t ~ t r - nx){ cos[2irf r (t - x R - nx)]} 
P 


where 


dsin0 


x = 


= differential time delay between 2 consecutive array elements 


2R 


x_ = — = round trip time to target; f = received frequency 
K c r 

For a static target, the received frequency f^ equals the transmitted 
frequency f^. However for a target moving at a radial rate the carrier 
frequency will be modified by the Doppler effect and 

f = f ± 2 ~ — = f + f . 
r c dt c c d 

where the sign depends on the sense of the radial movement of the target. 

Very high carrier frequencies are used. To process these, the received 
signals are heterodyned down in frequency. The signal received at the nth 
element of the array after heterodying then has the form 

v in (t) = p T (t - nx)cos[2Ti(f if + f d )t - 2r (f ^ + f d )nx] 


where the time origin is now taken at t = t , f , is the Doppler frequency and 

R a 

is the intermediate frequency of heterodyning. 
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If the maximum delay in the received signal across all of the detector 
elements is small in comparison with the pulse width, and if the doppler 
frequency shift is small in comparison with the carrier frequency (these 
are reasonable assumptions for real targets and existing arrays) that is if 


ni « T for all n. and 
P 


f, <K f 

d c 


then the approximate heterodyne signal at the nth array element can be written 
as 




v in (t) = p T (t)cos[2ir(f. !<: + f^)t - 2irna] 


if * d ' 


where 


and 


«- T «c + V “ T sln6 o 

c 


= free-space wavelength of transmitted carrier. 


only those 


will be completely processed where is the integration time for the proces- 
sor and Tp is the pulse width. 

With the electron gun scanning in a raster format, the input electronic 
signals can be arranged in the proper sequence to be compatible with the 
raster scan format by using N delay lines and an electronic switch. 


If the range of interest corresponds to a range delay 
signals., reflected by targets whose range delays satisfy 


- 


(T. - T ) 2R. 

__i E_ - i 


< T„ + 


C - R 


(T. - T ) 
i E_ 
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For now we will assume that the input data in the format of Fig. 3.2.1 
exists. Since only negative signal polarities can be formed by the electron 
beam, the most faithful optical transmission obtainable is 


+ 

2 

T(x,y) = l <5(y - n£)[l + yv (x)]P (x) 

N-l 
n= n 


Where P n (x) = - 0/2 

D 0, otherwxse. 


The distance x across the crystal is related to time t by 


x 


D % 

■zr t 


i 

is the integration time, £ is the vertical separation between scan lines 
in the input, D is the crystal aperture along the x direction, and y is 
the modulation index 0 <_ y <_ 1. T(x,y) is the ratio of transmitted/incident 
light amplitude. 

This image, and consequently the output data also, will remain on the 
crystal until it is electrically erased or until it decays through 
leakage currents (see Chapter 1) . Delay line modulators may use similar 
input data image formats for data processing, but the output is available only 
during the short time that the signals fill the delay medium. 


3.3 Electro-Optic Processor 

The optical system used is shown in Fig. 3.3.1. A choice must be made 
between the two possible light modulation schemes, amplitude modulation and 
phase modulation, described in Section 1.3. For reasons to be explained later, 
amplitude modulation of light is preferable. The electro-optic crystal trans- 
ducer is then placed between crossed polarizers. A general analysis of the 
signals and processing involved will now be presented. 



Mirror 



Fig. 3.3.1 Optical System for Radar Processing 
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The ideal optical processor illuminates the input (x,y) plane in Fig. 3.3.1 
with a collimated coherent light beam and passes the image formed through the 
integrating or Fourier transform lens in Fig. 3.3.1. The back focal plane of 
this lens contains the Fourier transform F(u,v) of T(x,y) 


F(u, v) = E 


where 


u = 


T( Xl y)e- 2,, J <ux + y y>dxdy 




(u,v) = back focal plane spatial frequency coordinates 

(x'jy 1 ) = back focal plane geometrical coordinates 

= light wavelength 

f = focal length of transform lens 

E = constant dependent on optical system 

F(u,v) = back focal plane light amplitude distribution 

With a crystal transmissivity ^ 

N-l 

2 

T(x,y) = l 6 (y - n£)[l + yv (x)]P (x) 

N-l in D 

n= t 


and an input IF modulated signal 

T 

0 

p 

where 


T. 


v in (x) = p T (j^)cos[2Tr(f if + f d )x ^ - 2 Tina] 


a - T(f c + £ d ) - fsine 

c 

Substituting T(x,y) into F(u,v) the output plane pattern is found to be 


oo N-l 


F(u,v) = E 


I S(y - ni)[l + yv ( x + ^dxdy 

N-l ln U 

n= r — 


rr . 
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substituting V^ n (x) F(u,v) becomes 


( f 


T T 

l 6 (y - n£)[l + P T (x ^-)cos{2n(f if + f d )x^-2ima} ]P n (x)e" 27tj (uidV y) dxdj 


D 


E 


I 6(- y-n£)P D (x)e- 27rj(uX + vy) dxdy 


( 1 ) 


+ E 


I <s(y - 


n£)Pj. (x) cos [2tt (f^ + f d )xj^ “2irna]e 2tt ^ ( ux + '^dxdy 
P 11 L 


where 


1 for |x| <_ X /2 


x - i 2 - d £ i> and p (x) , ' ; - p- 

P X 0 otherwise 

P 


The 1st term can be written as 

f CO, 


( 1 ) = 


l <5(y - n£)e 


-2Trj(ux + vy)P D (x) dxdy 


Cl 


6(y - n£)e" 2irjvy dy) 


„ , N -2irjux, 
P D (x)e J dx = 


N-l 


N-l 


SimrDu 


I e -2"jvnH , siajDu (1 + 2 J cos 


ttDu ~ ttDu , 

N-l 1 


simrDu 


ttDu 


2 - n 

N-l 

f |/C T . -A ”1 sin C 2 )ttv£ 

!i + 2coJrf J y* v M— — - 

l-* -1 simrv£ 


simrDu ^ + sinNrrvl. _ simrvg. ^ _ sinuDu sinNiTv£ 


irDu 


simrv£ sinrrv£ irDu simrv£ 


while the second term can similarly be manipulated into the form 


I <5(y “ n£)P x (x)cos[2Tr(f i<f + f^)^ x - 2irna]e 27r j Cux + vy) dxdy 


’if d D 


( 2 ) 


( 2 ) = 
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1 

2 




I My - 


n£)P x (x) [e 
P 


2 TTj ( f if + x “ 2lT 3 na -2-rrj (f ±f + f d ) 


+ e 


D X+ 2 ’ Jn “]e- 2,t;i<ux + Vy) dxdy = 


-hi 


S(y - n£) e - 2 ,I 3na e -2ii jvy^ 


P Y (x)e 


2 *3«if + Vd 1 V 2 "J ux dxl + 


+h i 


r 00 


T. 


6(y - „H)e 2,, j n V 2,,Jvy dy 


-2lTj(f. J . + f,)— X 0 . 

P v Me l£ d D e- 2,,: > ux d x ] - 


N-l 


N-l 


_ sinirX (u - - — ) 2 

1 r £ 


“f i 


K , ' | e -2irjn(a + vZ) + slnTrX p^ u + ^ e -2irjn(a - v£) -j 

nX (u - T- ) N-l uX (u + r- N-l 

p A ~ ~ v \ ~ — 


, sinTrX (u - -r— ) 

J- r P Kr 


•i i 


ttX (u + — ) 

P *x 




+ 2cos — - ir(a + vZ) 


sin( ^ ^ ) tt (ct + vZ) 


sinxCa + v£) 


-} + 


sinirX (u + — ) m 4. i 

E -r— x~ {i + 2cos - — 7 , — Ti(a + vJl) 

ttX(u + j-) 

P X x 


sin( N ^ )iT(a - vZ ) 


sinn(a - vZ) 


-}] = 


, sinirX (u + — — ) . „ , . sinirX (u + ■— ) . 

„ 1 r p 1, sinNir(a 4- vZ) + p Ay sinNuCa - vZ) , 

2 v , lx sinn(a + vJl) v , . 1 ■. simr(a - vZ) J 

X (u - — ) ttX (u + t— ) 

P X v P 


where 


X X = + f d ) • « ' sl " 6 o 


Combining these expressions, we find 

F(u,v) - (1) + (2) - E SiaiSH. siaSjvl + 
* jr Du sinrrvS. 


„ sinuX (u - - — ) sinNir(a + vJl) 

+ f E— Jx 

irX (u - - — ) siniT(a + vZ) 

P *x 


(*) 
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+ f 


simrX (u + - — ) 

n A 


-x- 


nX (u + — ) 
p Ax 


sinNirCa - v&) 
simr(a - vSL) 


Fig. 3.3.2 shows this output pattern. 

With the variable change vH = - sin 0, we obtain 

c 

. _ sin(N7T^— sin6) 

F( u ,e) = E 


ttDu 


sin(TTr— sin0) 
c 


+ 2 


c, sirnrX (u - 7— ) sin[Nrv— (sinQ - sin0 )] 
E P k k ° 


ttX (u - - — ) sm [ 7T- — (sm8 - sin0 )] 
p A A o 

r x c 

sinirX_ (u + 7—) sin[NT!Y— (sin0 + sin0 )] 


_k_ 


2 nx (u + f-) 
P k 


sin(iT— (sin0 + sin0 ) ] 


where the dependence on 0^ has been emphasized. 

In this derivation, we have assumed that T. >> T or D » X which is 

1 P P 

the case for real time processing of radar signals at horizontal TV scan 
rates. The first term in the F(u,v) expression is the zero-order 
pattern resulting from the constant in the 1 + v form of T(x,y). This 
term yields no useful information. 

The second term (*) is the positive first order fringe; the position 
of its lowest order maximum is at 


“i - r - r <£ if + V 

x 


_ cx - d 
V 1 = “ 4 ~ AX 


sin0 


c 

As shown in Fig. 3.3. 2,f^ is assumed to be greater than in absolute 
value so that there is no indetermination i n deciding which first order 
spot to choose, the one corresponding to the positive u axis is chosen 




f f 






Fig. 3.3.2 Output Transform Plane Pattern for the Optical Processor 
of Fig. 3.3.1 with the Input Format of Fig. 3.2.1 
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From these expressions and Fig. 3.3. 2, we see that we can obtain 
both range rate and angle information about the target from the position 
of this output plane maximum. The linearity of the system allows multi- 
target processing capability also. 

Furthermore, the amplitude profile of this first order maximum will 
be the same as the amplitude profile of the antenna beam itself. (Ref 9) 

2 sin^[irN y-(sin0 - sin0 Q )] 

| F(u = — , 0) 1 « i V(9) | 2 = K(0) r -~ 

x sin [tit — (sin© - sin0 ) ] 

A O 

C 

where K(0) represents the radiation pattern of an individual element in the 
arrayJ.K(0) usually varies only slightly from unity for the values of 0 
of interest, as shown in Chapter 2. This proves that the output of the 
optical processor corresponds to the output of the ideal processor defined in 
Chapter 2. The Doppler resolution obtained is the same as that of an ideal 
processor with integration time T^. 

The result just obtained, indicating that the angular and frequency re- 
solutions of the processed signal are the. same as the ideal output resolutions, 
depends on the implicit assumptions that the optical system is diffraction 
limited and that the pattern on the input crystal is a faithful representation 
of the idealized IF radar signal. For the second assumption to be valid, the 
crystal- resolution must be sufficient to reproduce the T^(f^ + f^) cycles of 
signal along each line in the x direction in proper relative phase to the sig- 
nals along the other lines. For processing of an N element linear array, this 
condition requires a crystal resolution of at least N lines in the y direction. 
These minimal values must be weighted against the attainable resolution of the 
target crystal. 

Pockels tubes operated at the Curie temperature exhibit resolutions of 1000 


lines in 5 cm square crystals. This means that intermediate frequencies up 
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to 


1000/2 

50xl0~ 6 


= 10 MHz 


can be faithfully represented on the crystal. The, restr iction in N allows 
us to process returns for arrays of up to 500 elements, by far larger than 
any existing linear array. 

Another remark about T(x,y) must be made. The delta functions 6(y-n£) 
account for the fact that the time history from each detector element is writ 
ten at discrete vertical locations y = n£ on the crystal, where n = 0...N and 
SL is the separation between vertical channels on the input crystal. This 
approach simplifies the summations and computations involved. 

However, if we consider the finite size of the spot (1 mil) and assume 
the cross-section transmission of a horizontal line to be given by 

g(y) = p (y) =/ 1,0 - |y ' -f 

y s \j) otherwise 


where Is is the spot size, the value of F(u,v) formerly derived must be modi- 
fied to 


„ sin nDu „ , N sin N ttv£ 
F(u,v) = E G (v) 


ttDu 


sin Ttvi 


♦I 


sin irX (u - — ) 
P 4 


TtX <u - 
P X x 


x n f \ sin N i(a + v£) 
1 ^ sin it (a + v£) 


+ f 


sin irX (u + t — ) 
P A 


*X (u + —■) 
P A 
r x 


x r /\ sin N TT(a - v£) 
1 v sin it (a - v£) 


/ CO 

g (y ) e ^ 2TrVy dy is the Fourier transform of g(y). 

•OO 
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A numerical example will be most instructive here. 
G(v) is given by 

AG 


The 3 dB width of 


1 

'3dB = s 


25 x 10 


-6 


= 4 x 10^ m ^ 


while the beamforming area for 6 = ±90°, considering — = 0.5, is given by 


v 

max 



2 x 10 


-4 


5 x 10 3 


-1 

m 


from which we find that G(v) is almost constant along the beam forming area for 


6. The approximation of g(y) by 6(y) is then valid. 

3.4 Data Format Considerations 

Ttoo methods have been considered for producing the format depicted 
in Fig.. 3.2.1 and we now devote some attention to the advantages 
of each scheme. 

3.4.1 Time-Delay Multiplexing 

In this implementation, the time stories of the range gated IF signals 
for every element are delayed an amount of time which is a multiple of the 
horizontal scan time T^ of the electron gun prior to being written on the 
crystal sequentially element by element. Fig. 3.4.1 shows this situation. 

This method requires N delay lines , an electronic switch and support 
electronic equipment. The pulse repetition rate of the transmitting 
antenna has to be lower or equal to the 30 frames per second video rate. 

3.4.2 Time Sampling 

Another implementation consists of sampling the IF signal present across 
every array element during the total video frame time. Fig. 3.4.2 and 3.4.3 
illustrate this input format. This procedure is similar to time multiplexing 
but in this case the electronic switch will sample each of the N elements 
of the detector array at the rate of 1 sample per element every 50 ps. 

The electron gun now scans vertically at the normal horizontal scan rate 


of TV. 


-4 


! 

i * 
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. Fig. 3.4.3 shows the actual detector signal and the result after 
sampling. For the sampled signal to accurately represent' the real sig- 
nal, the frequency and pulse width of the signal are severely constrained. 

14 

For TV writing rates, f = -z ^ = 10 Hz is the maximum frequency 

max / x du y s 

which can be accurately represented on the crystal, this corresponds to 
two samples per cycle. The pulse width transmitted by the antenna must 
now be larger than 100 ys , if a one cycle cosinusoidal signal at the maxi- 
mum frequency is to be written. 

It must be remembered that now the electron gun is scanning the 
crystal vertically , at the rate of 1 vertical line every 50 ys. This 
scheme again requires the least support electronics. Its disadvantage in 
radar systems is the large pulse width required. This scheme can be 
shown to be quite applicable and practical with sonar data however. 

3.5 Amplitude versus Phase- Modulation of Light 

.. Amplitude modulation and phase modulation of light present both 
advantages and disadvantages. Phase modulation of light presents a 
higher diffraction efficiency than amplitude modulation; furthermore, 
the phase modulation is linearly related to the voltage applied across the 
crystal (as shown in Section 1.3.2). While the amplitude transmittance 
varies as a sinusoidal function of such voltage (Section 1.3.1), being 
approximately linear for 5 <_ 30° with error less than 10%. 

Phase modulation of light requires a larger optical system bandwidth 
than amplitude modulation, for the same amount of distortion. Furthermore, 
it introduces greater noise consequently reducing the signal-to-noise ratio 
at the output plane, and it is more sensitive to irregularities in the 
crystal and windows surfaces. Based on these facts and observations performed 
for both modulation types using the same input patterns, the choice to use 


m?” 
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amplitude modulation for this particular kind of processing was made. 
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CHAPTER 4 

ELECTROOPTIC PROCESSING OF TWO-DIMENSIONAL PHASED ARRAY ANTENNA SIGNALS 

4.1 Introduction 

Two dimensional phased array signal processing consists essentially 
of correlating the phase slope corresponding to the planar waveform im- 
pinging on the array proceeding from a reflecting point target 
at a range R at a direction given by two angles refered to a coordinate 
system consisting of the array itself and the boresight direction. All 
that was said about the linear array applies here, with the added 
requirement that the output of the ideal processor must now yield information 
on two angles at a time. Two input formats will be considered, their 
requirements compared and their on-line implementation discussed. 

» 

4.2 Full Array Processing 

We consider the receiving array depicted in Fig. 4.2.1 where the 
received waveform has a plane equiphase- front whose normal is defined by 
the two angles <^. An array consisting of a uniform rectangular grid 

of M elements spaced at a distance AX in the X direction and N elements 
spaced at a distance AY in the Y direction. The input format used consists 
of dividing each of N lines in the electrooptic crystal into M sections 
or cells. The mth cell of the nth line contains the representation of 
the signal received at element (m,n) of the array after amplifying and 
heterodyning during the integration time T^. 

Let this i.f. signal be represented as 

v i (t) “ Pt (t - V cos[2 " (£ if + V' - 2 "< f c + V T „ m ] 

nm P 
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where 



differential delay per element 
d 


nm 


c 

distance from antenna element (n,m) to the impinging plane 
RF carrier freauencv wavefront 


IF frequency 

Doppler frequency shift (f^ <« f and f^ « f^) 


When r << T for all m, n which is a reasonable approximation in most 
nm p 

cases we have. 


where 


v inm (t) “ P Ip <t)cos[2„<f 1£ + f,)t - 2™(f c + 


d mAX cos <J>„ + nAY cos <j> 
nm _ X Y 


mn 


follows from the geometry of Fig. 4.2.1, and 

mAX nAY 


(f + f ,)t = f t 

c d mn c mn 


T~ cos ♦x + T 


cos 4 > , 


mAX . nAY . 

= sxn0 Y + sinB 

A A. A I 

C C 


<f>. 


x and <f> are related by 

y e v = 90° - 4> 

x X 


0 = 9 0° - <f>Y 


if we define two a parameters 


AX 4 A 

“x ■ — sln0 x 

C 

AY 4 A 

“y = r sine v 

c 


we can rewrite v^ as 
nm 


v i (t) = P T (t) COS f 2 7T ( f . £ + f d )t 
nm p , 


2Ti(ma x + na y )] 
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With the input format depicted in Fig. 4.2.2 the optical transmission 


function is given by 


N-l 

2 

T (x,y) = l 6(y - y n )P D (x)[l + p 


n=- 


N-l 


M-l 

2 

l 


x - x 


m=- 


M-l 


P D/M (X 


x n )v i 


n 


nm 


D/M T i^ ■* 


where X and Y are the geometrical coordinates of the array (Fig. 4.2.1) and 

x and y are the geometrical coordinates in the input crystal plane (Fig. 4.2.2) 

y and x are related to X , Y by two scaling factors k and k 
n m mm mt x y 

Xmn mAX D 

x =*—. = — — = m — 

m k k M 

x x 

y = Ymn = mAY = n £ 

n k k 

y y „ 

As shown in Fig. 4.2.2, M is the number of channels, l is the vertical spacing 
between them P 1 , (x) = 1, ]x| .<_ *^2M and 


M 


0, otherwise 


D is the crystal dimension across the x axis. Using the optical 
system shown in Fig. 3.3.1 a two dimensional Fourier transform of this 
input plane amplitude distribution is taken. 

The output plane pattern is given by 

f 00 I 


F(u,v) = E 


= E 


r 


T( X ,y)e- 2 ^ (ux + v y ) dxdy 


r ” 


N-l 


M-l 


n=- 


^-l 6(y ' n<1)P D <x)(1 + " I M1 P D/M (X - m 


m=- 


T M xT.M 

P X U “D" " “T i )c°s{2^[(f if + ~ oT ) - (ma + na ))} 

P A. 1 


e -2jj(ux + vy) dxdy 


where 


_ T 

A p M T ± 
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With T, < T , the P„ pulse is absorbed into the P , . 
i — p X u/ci 

» n 


Thus we obtain, 


decomposing the cosine term into the sum of two imaginary terms, 

N-l 

— 2irj (ux + vy) 

dxdy 


F(u,v) = E 


r 


l 6(y - n£)e 
N-l 


N-l 


+ yE { 


— 2iri i 


2 -'TinoY 

i -2ujvy £ 6 (y — n£)e dy} 

N-l 


r m 


_2„ 1x • T i M 

P D ( X )e 2 ”^[u - (£ lf + £ d ) -i- )} 


, Yp , „ X+ T i( f if +f d )) 

( ,A p d/m (x - m 5 )e > 


M-l 

2 

r 00 


N-l 


+ yE { 


e -2irjvy' l Hy _ ni)e 2^ dy} . 
N-l 


P D ( X )e- 2 "J x t u + (£. f f f > ^ ]} 


M-l 

{ Jl W* 
2 


N-l 

m I <5(y " ni)e 2 ^ na Ydy} 


N-l 


Expression which reduces to 


x T . sinNTrJlv sinirDu , 

F(u,v) = K : — + 

o sinirilv uDu 


+ K 


simrN(a y + JLv) simr||(u - ) sinMTi[a_ + ^u] 


x K 


1 simi(a + i v) D. 1 . . r , D , 

x "m (u - S1I "[“ X + if 1 
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+ K 


simrN(a - &v) simt|(u + ) sinMir[a - ^ u] 

X — r : ^ 


1 sln,< € ‘ tV) 4(u + ^> In^-S-l 


where 


* =J5 i 

x KT. (f. f + f d ) 


The second term in F(u,v) corresponds to the first order positive fringe 
of interest. This distribution is shown in Fig. 4.2.3 

It can be seen that this term corresponds to the antenna pattern in 
the v direction, with grating lobes occuring at 


v 1 = - y- + j, n = 0 , +1 , +2 , . . . 

. D, i . 

siniT M (u “ r 0 

In the u direction, two terms in u appear. The term — corresponds 


D, 

*M (u 


-) 


to the Fourier transform of the input cell in the x direction, is non- 

repetitive, presents a maximum at u = and varies smoothly over its 

v D 

half power width which is given by Au = The term sinMu[a 4- ^u] 


JL 


M 


sinufa + rju] 
X M 


corresponds to the antenna pattern on the x direction, is repetitive and 
presents main grating lobes spaced at values, of u given by 

u = - — a + — m, m = 0, +1, +2,... 
m D x D — — 


It can be seen that there are N resolvable beams in the v direction and 
M in the u direction. The angular information corresponding to 0^ can be 
extracted from the position of the zero order grating lobe in the v 
direction 

fSr 

v i " & 

Vc 

thus sin 9,, = = 

x AY 


jr v i 




i Fig. 4.2.4 Cross-section of the Output Pattern across v*=v^ showing the cell 

i pattern and the sampling term. 

1 - — / 
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To obtain the information corresponding to 6 we observe that the zero 

A 

order grating lobe of the sampling term in the u direction should fall 
inside the 3 dB width of the term corresponding to the input cell, for 
all values of 0^ of interest. This can be achieved by selecting f if T i 
to be an integer n^. The position of the n^th order grating lobe of the 
sampling term now occurs at 


u 


1 


(n. - a ) 
1 X 


M 

D 


and the maximum of the cell term occurs at 


1 M . M 

U 2 A D f if T i n i D 

x 


Figure 4.2.4 shows these peaks in an output plane cross-section 

The angular data corresponding to 0 is thus contained in the position of 

A. 

n ± th order grating lobe of the sampling term 

U 1 53 ^ f if T i ~ ? x- D 


from which 


a A £(u. - f T . ) A 

„ y c M 1 if i c 

sin0 v = — ? — 

X Ay Ax 


It has been shown that the optical processor with the input format de- 
scribed functions as the ideal processor for two-dimensional phased array, 
presenting a resolution of M cells in the u direction (corresponding to 
0^) and N cells in the v direction (corresponding to 0 ) , in accordance 
with the Rayleigh resolution criterion. However the input plane requirements 
are very severe, independently of real time operation for the full array 
format to be implemented with the current crystal resolution. Indeed, 

N - 100 lines across the vertical dimension are easily accomodated; in 

the horizontal direction, however we must accommodate M - 100 time 
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histories. For f. - = 5 MHz and T = 1 ysec, this corresponds roughly 
if p 

to the full crystal resolution, using only two signal samples per cycle 
and we recall that we must also maintain the relative phase of each cell 
accurate to one part in N. 

Real time processing of the full array would always require -at least 
M x N delay lines, amplifiers and heterodyners . 

Even with the low degree of feasibility of this format, its study is 
of interest because modifications of it will lead to a far more feasible 
input format. 

4.3 Thinned Array Processing 

4.3.1 Input Format 

The input resolution requirements for a full array processing can 
be relaxed by processing only the returns corresponding to 

some of the array elements, that is, by taking a thinned array. Again 
in this case the time histories of various detector elements are recorded 
in separate channels along the crystal. However, only one element in 
each row and column of the receiving array is considered, i.e., a total 
of N elements assuming a square planar array of N x N elements. Fig. 

4.3.1 shows the thinned array geommetry. 

In the input plane format, the y coordinate of each channel, y^, is 
directly proportional to the y| geometrical displacement of the ith sensor 
from the center of the array. Similarly, the center of the ith received 
waveform is displaced by an amount x^ from the optic axis corresponding 

or 

to the displacement of the sensor from the center of the array. Figure 

4.3.2 shows this input format. With the appropriate scaling, time and 

x displacements are comparable. If we represent a length D of. crystal as a 


'W 
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Fig. 4.3.1 Sampling and Notation for the Thinned Array Input Format 



channels 


Fig. 4.3.2 Input Format for 


a Thinned Array Planar Phased Array Processor 
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time T, time and displacement in the input plane are compatible and for the 
spatial signal the relations are 


T D 

t = x D ; X= T fc 


To avoid retaining additional unnecessary terms in the ' computation, the 
complex signal (t) corresponding to the real signal (t) is considered. 
This will result in fewer terms in the output plane which correspond 
exactly to the terms of interest to determine angular information. The 
complex form of the received signal at element i of the receiver array 

i 

prior to heterodyning is 


V ± (t) = s T (t - T i )exp[2iTj (f c + f d )t - 2mj(f c + f d )T ± ] 

where s_ (t) is the signal envelope, assumed of temporal duration T 
T P P 

(Requirements on s(t) will be considered later. After heterodyning 


and assuming << T , 


V i ^ = S T <Oexp[2ir;j(f + f d )t - 2irj(f c + f^xj 

P 1 

The spatial signal is given by 

V ± (x) = s x (x)exp[2Trj (f if + f d ) ^ x - 2irj(f c + f^x.^] 

P 

where X = ^ T . Figure 4.3.2 shows this input format. The time delay x. 

pip r i 

for element i, of geommetrical coordinates X d 8 ,Y. 8 iri the real array 


d X 8 cos + Y. g cos<fc . 

__ i _ i X i y Y _ 

* _ " " *“ 


i c 


k sin0 

8 -5 JL + v S 

c J ± 


k sin0 


where ^ , •J’y » 0^ and 0^ were defined in paragraph 4.2 and 


1ST- 



- 61 - 


X 8 Y 8 

k = — =— and k = are scaling factors which relate the geometries of 

x x 8 y v 8 

x i y i 

the receiving array and the input plane. can thus be written . as 

k sin0 k sin0 

V ± (x) = s^(x)exp[2TTjf^ f x - 2irjf ^(x^ 8 X — + y ± 8 V g— ~ )] 


where we have assumed 


and where 


f » f . , » f 

c if d 


£ h - f if 5 “ < £ if + V F 


p p 

After being delayed by an amount x^ corresponding to its coordinate 
on the array such that 

X 8 

„ 8 L_ 

X i " k 

x 

the form of the signal written on the crystal for the ith element is 


V i (x) = s^ (x - x i 6 )exp [2iTjf^ f (x - x^) - 2ir j (a x. 8 + a y. 8 )] 
P 1 1 x l y i 


where 


f k sin0 k sin0 

a = -. c x X = — X 

x c A 

c 


f k sin0„ k sin0„ 
a = c y Y = Y 

y c A 

c 


The optical transmission function will again be of the form 


N-l 

2 

T(x,y) = E l [1 + pV (x)]6(y - y 8 ) 

Mzi 

. ■*■*“"* 


2 
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The Fourier transform of this amplitude distribution is taken by the 
optical system described in Fig. 3.3.1 . 

The first order fringe of interest in the output plane is given by 

N-l 


F 1 (u,v) = E 


2 2u j [ f ! _ (x - x. 8 ) - a x. 8 - a y. 8 ] 

g N J if i xi y J i 1 


-CO 1 =- 


l <s(y - y, 8 )e 

N-l 


e -2nj(ux + vy) dj , dy 


N-l 

2 t 

-yE l i 

N-l 


_ . -2irjy. 8 a 

6(y - y i g )e _2 * ;|Vy e 1 y ) 


■ f 
{ 


2ifj[f! i: (x - x. 8 ) - a x. 8 ] 0 . 

„ / gv if x x x J -2ttjxu 

s (x - x. )e e J 

A 1 

P 


dx} = 


N-l 

2 t 

yE l { 
N-l 


- -2*Jv<y - y 8 ) -2xjy 8 (v + a ) 

5(y-y 8 )e 1 e ? 


dy} 


s x (x. 


x. 6 )e 2,T;i[f i£ <x - x i S) 

X 


.- (a + u)x. ] -2irju(x - x. 8 ) 


1 dx} 


N-l 


2 -2irjy 8 (v + a ) 


y }{S (u - f» f )e 


-2irj (u + a )x. 8 


x x 


^ I {e 
_ NKL 
2 

N-l 

2 -2Trjy 8 (v + a ) -2irj ( w + f' + a )x. 8 

PE l {e 1 y }S(w) {e lf x 1 } 

N-l 
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N-l 

2 -2Trjy 8 (v + a ) -2irj x 8 (w+ f! + a ) 
l e y e 11 x 
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where u and v were defined in Chapter 2 , and w is given by 


W " U " f If 


and 


S(w) = 


, . -j2TTXW, 
s x (x) e J dx 


is the Fourier transform of the transmitted signal envelope s(t). 

T P 

The grating lobes of F^(u,v) occur at points (u,v) in the output 
plane which satisfy 


2tt{ 


[y 8 (v + a )] + [x 8 (w + f* + a )]} = m2n, 
i y 1 it x 


where m is an integer. The desired angular information (0 and 0 V ) can be obtained 

X. I 

by considering those points in the (u,v) plane for which 


y^ (v + a y ) = n/, an integer 


p * X 

x ± 6 (w + a x + = nu , an integer 


± = _ Nzi o 
x 2 > • • • • 2 


By choosing f!^x^ 8 to be an integer for all i, these equations reduce 

to 


y i 8 (v + ct y ) = n ± y 


x i 8<u + “x> - n i X 


for all i. 


The zero order grating lobe is obtained for n^ y 


= =0, yielding the 


coordinates 


k sin0 v 
y i 

v = -a ■*— r 

1 y A 


k sin0 
x Y 

u- = -a — 

l x X 
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thus 


sin 6^ 


c 

r v i 


X 

sinQ x = - IT U 1 


Two requirements for this processing to' be valid are: 

1) The half power width of the transmitted signal spectrum S(w) must 
be large enough to accomodate all possible beams for 0 within 
the beam forming cone. 

2) It must be proved that grating lobes of F^(w,v) ’ occur only at the 
points for which 


y i 8(v + V = n i Y 



(u +- a ) = n. X 
x x 


for all i. 


The design parameters must be thus chosen so that both of these conditions 
are satisfied. It is also a design objective to keep the sidelobe level 
of F^(u,v) as near as possible to the sidelobe level corresponding to the 
full array aperture. 

The crystal resolution must now be able to accomodate only (f_^ + f^)T^ 
cycles of IF signal across the x direction, which is a factor of N improve- 
ment over the full array requirement. 

4.3.2 Choice of Design Parameters 

The following parameters are considered 


Carrier frequency f = 5 GHz 

-2 

Carrier wavelength X £ = 6 x 10 m 

p p —2 

Spacing across array elements AX.^ 6 = A.y^ = L = 3 x 10 m. 
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Scaling factors = k x = 100 
Intermediate frequency f^ = 3.333 MHz 

Transmitted envelope s(t) = 13-bit Barker coded signal with = 13 ysec. 

Number of array elements = 4900; N = 70 

Cone of beam forming 0 <_ 30° 

X 

e Y < 3°° 

The following points will prove the feasibility of the processing: 

a) f is large enough so that the quantizing of x^ in units of 1/f!^ 
leads to negligible errors. Indeed: 


f! = f ~ = 3.3 x 10 6 x 50 — - 10 — - 
lf lfD 5 x 10~ 2 


= 3.3 x 10' 


1/f ! £ = 3 x 10 
if 


-4 


x 6 , 3 x IQ' 2 _ -4 

i min 100 J X XU 


thus the quantization is exact even for the smallest x^ . 

b) The maximum value of x^® is much less than the aperture width of the 
crystal so that nearly complete waveform stories can be recorded 


there without significant loss from edge effects 

... . . --' - 2 

Indeed, x. 


70 x 3 x 10' 


l max 


100 


= 2.1 x 10 2 in. 


With 


-2 

v _ 5 x 10 ,„-6 , „ , rt -2 

X - — T t- 13 x 10 = 1.3 x 10 m. 

p T i p 5 x 10" 6 


_2 

complete time histories can be recorded on the 5 cm (5 x 10 m) 
crystal aperture. 

c) The angular resolution desired in beam forming and the scale factor 

Or or 

relating x^ and y^ & to the sensor locations in the array must be 
chosen so that adjacent beams are not more closely spaced in the 


Fourier plane than the diffraction limit allows. The diffraction 
limit in the system is approximately 


1 

D 


1 

5 x 10' 


20 m 


The attainable resolution in beam forming corresponds to the half- 
power beam width of the main lobe in F^(u,v) which is 


r\ 

NAx^ 70 x 3 x 10 /100 


100 . , -1 
23 “ 47 ” 


The required resolution in a^, a thus exceeds the diffraction limit 
of the optical system. 

d) The spectrum of the transmitted signal must be broad enough to ac- 
comodate all beam positions within the beam forming cone for 0 30°. 

For 0 = 30°, the value of a is 

X X 


„ _ 100 x 1/2 10 -1 

“X rcax - , ,„-2 " 12 ~ 830 m 

6 x 10 


The choice of signal s(t) must take this fact into account as well as 
simplicity of implementation when simulations in real time must be 
; run. Such conditions are met by envelope coded signals such as Barker 
coded signals. Plots of the spectrum S(w) for 7, 11 and 13 bit 
Barker coded signals are shown and the choice was made of the 13 bit 
coded envelope with a duration of 13 ysec which affords the smoothest 
variation in the area |w| <^0.8 x 10 ^ Hz corresponding to 800 m ^ . 

Fig. 4.3.3 shows such plots. 

e) The choice of sensor distribution must be done in such a way that 

ambiguous beams do not appear in the output plane for values of (u,v) 
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other than those considered in 4.3.1. We shall now prove that such 
is the case when the choice of sensors over the array( one for each 
row and column) is randomly performed. 

The equations to consider are 


x i 8 (” + \> ■ n i ! 


y ± (v + a y ) = V 


For a value n^ 4 0, for instance, we have 


x. K (w + a ) = n. 
i x l 


this is 

but since 
we have 


x.^w+ x.^a = n . 5 
i i x x 


x.^f!, = n., for all i, 
x if x’ 


n. n. 

x , x x 

-rr~ w + — a = n. 

fl* f' x x 

xf xf 


fT“ ( w + a J = n. X 
if x 1 


n 


yielding w = -a^ + — — f!^ for the possible maxima of F^(w,v). Two cases 

i 

now arise: 

I) The ratio of n^ and n^ depends on i. In this case, the minimum 
variability in the value of w is given when n^ =1 and 


n. = n . = x. x f' = 

i x max x max if 


i.05 „ „ . n 3 

— x 3.3 x 10 

10 ^ 


= 35 
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thus Aw 


mm 


= 33 x 3 . 


3 x 10 = 


100 

1.05 


= 95 m 


-1 


Thus even the smallest variability possible for Aw exceeds the output 
plane resolution of 20 m ^ and in-phase contributions from all or 
even several of the sensors is not possible within a single resolution 
interval in the Fourier plane. 

II) The ratio of n^ and n^ is constant and independent of i. In this 

Cr 

case, since x. f' = n. , we have 
’ i if i 


n 


(w + a ) = n X 
r if X 1 


which implies 


w + a 

— “ n , an integer 

if 


and thus 


w = -a +nf', n=+l, + 2 .. 
x if’ — — 


are the positions of the grating lobes other than the zero order 

lobe (similar reasoning can be applied along the v direction) . 

^ 3-1 

With f^ = 3.3 x 10 m , no ambiguous grating lobes occur within 
the beam forming angle of 1 9 1 _< 30°. 

It remains to be proved that no maxima of F^(w,v) occur for these 

V 

directions — for which 
w 

y. g (v + a) + x g (w + a ) = n , an integer. 

• - y i x ■ l 

. In such a case, we have 


y^ 8 (v + a ) = n ± - x i S (w + a x ) 


n. 


g 


v + a • = 

y 


— (w + a ) 

g x 
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this is 


-it! + it--ii! 

V = " a y y 8 “ x y 8 y 8 

y i y i y i 


w. 


Three cases arise now: 


I) n ± - 0 


Then the set of straight lines 

. n i 

v = -a + 

y „ 8 


— (w + a ) 
g x 


reduces to 


-a - 

y 


— (w + a ) 
g x 7 


x. c 

1 


which, for randomly chosen sensors across the array, this is, = f(i) 


interset only at the point 


v 1 = -a 

1 y 


w, = -a 
1 x 


thus in phase contributions for all of the array elements are possible 


II) 


only on the zero order lobe. 

n -. 

= R, constant 


g 


Then we have 


g 


x, 8 

v = -a + R — (w + a ) 

y v e x 

y i 


which, for — — = f(i), represents a set of straight lines intersecting 
at the point 

w = -a 


v = -a + R 

y 


It remains to be proved that R can take on only the values + m y- 

X/ 

where l = |f- and m = +1, +2... 

y • ' 
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and since y ® is quantized in units of .£, it turns out that the pro- 
duct y 8 ^ is an integer. 7 - plays here the role of f* when ambiguities 

lx, H 

across the u direction were considered and thus this problem reduces 


to one already solved. 


n. 


III) = R. , some function of i. 

g i 

nr w 


In this case we have 


x. 


g 


v = -a + R. - — (a + w) 
y i g x 

i 

which is a set of straight lines that cannot intersect at a single 
point on the (u,v) plane, what would require R^ to be constant and 
independent of i in the case of sensors chosen randomly across the 
array. Thus in phase contributions from all or several sensors is 
not possible outside the grating lobes formerly considered and the 


proof is complete. As a conclusion, Fig. 4.3.4 shows the position of 

the grating lobes for 0 = 0 = 0 °. 

X . Y 

4.3.3 Choice of Array Sensors 

The criteria to be used for selecting the location of the elements 
in the thinned array are 

1. The resolution in the transform pattern of the entire N x N array 
should be maintained 

2. The sidelobe levels in the transform pattern should be as low as 

possible. 

2 

When N elements of the array were used, the optical beam could 
2 

be steered into N resolvable positions. For the thinned array of N 
elements, in order to avoid the observation of false beam positions, 
it is necessary to minimize the peak power in lobes other than the main 
grating lobes. This problem has been the object of much research and up 
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to date no better approach to solve the problem is known than a pseudo- 
random choice of the elements across the array (as proved by intensive 
computer research of this specific thinned array pattern by A.P.L, 
Reference 8 ) . Such pseudorandom choices have been proven 

to yield sidelobe levels 10 dB down from the main lobe peak power 
against the 13 dB attainable for the full array beam. (Ref. 8, 10, 11) 
A computer program has been written that automatically selects 
the elements in a planar array in a pseudorandom mode so that several 
configurations can be evaluated and the best pattern orientation chosen. 


• . 
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CHAPTER 5 . ... 

. r 

COMPUTER INTERFACE *• . 

"* • 

.* 5.0 Introd unction . ‘ 

• • 

One of the major accomplishments. in the project is. the comoletion of the 
computer interface. This provides the connecting l.ink from the optical 
processor to a digital computer. The resultant hybrid processor was first 
proposed in April 1972 at the Optical Computing Symposium in Darien, Conn. 
It extends the versatility and processing power of the system further 

beyond those of other optical computers. It enables one to utilize the 

. » 

optical processor for those operations for which it is best suited and 
the digital computer for those areas in which it performs best. 

The interface is constructed from modular computer functional blocks, 
known as Register Transfer Modules, RTMs that were designed here at CMU. 
Reference 12 provides some of the design details of these 

computer modules. The forthcoming Sept. 1973 issue of IEEE trans.on 
Computers contains'' a brief description of this interface, Ref. 6. .A 

complete description of the work is one; of the many papers that must be 
writte*n this summer and submitted for publication. 

The modules used were loaned from ARPA funds and should at some 
point be purchased. The PDP-11 to which the system is interfaced is also 
shared with the RTM program in the department. If free money ever appears, 
ve should purchase our own since its workload increases monthly and we 
could move much faster with our own unit. 
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5..1 OUTPUT INTERFACE 

An extensive discussion of this optical-digital interface is in- 
cluded, as it appears to offer considerable interest. The interface 
essentially digitizes each TV frame into a number of cells or picture 
elements, 4 x 4 to 128 x 128 cells at a rate of 33 frames/sec (standard 
TV rates) . The intensity of each cell will be either a 1 or a 0. This 
decision being determined by whether the light intensity within each 
cell exceeds some progr ama ble threshold level. The resolution of each 
frame is under computer control. The prime use of the system is in 
locating points of light in the output correlation plane of a real time 
optical processor, this search being done automatically without the 
need for human intervention. Figure 5.1.1. illustrates sample results 
of this digitization process. 

The block diagram of Fig. 5.1.2 shows the system with the interface 
enclosed by the dashed lines. A vidicon will be continually scanning 
the output correlation plane, and its composite output signal fed to a 
sync separator. For timing reasons, scanning the output correlation plane 
and its vidicon. is driven by an external sync generator central to the 
entire optical data processing system. The TV signal then enters the 
video thresholder which digitizes each horizontal scan line, according 
to programmed horizontal resolution information. Timing for the threshold 
er is provided by horizontal and vertical blanking pulses from the sync 
separator. v 

This data is then buffered and transferred to the digitized frame 
generator which, according to vertical resolution information, combines 
data from successive horizontal lines to form the final picture cells. 

The contents of these picture cells are then transferred to a PDP-11 

which can analyze and display the data, and, conceivably control the 
entire optical process. 


6 - 







.Fig. 5.1.1 Illustration of TV frame digitizing 
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5 « 2 VIDEO THRESKOLDER 

A block diagram of the thresholder is shown in Fig. 5 / 2. 1 A. portion 

of the horizontal scan line (t seconds) is examined. This corresponds 

c 

to one cell size. If the voltage level of the signal (corresponding to 
light intensity) exceeds the threshold v^, during the time, a 1 is placed 
in the shift register shown. After 16 x c portions of the line have been 
examined, or the end of the horizontal line has been reached one 16-bit 
data word is transferred to the digitized frame generator. 

The time x c is determined by the value of the frequency control 
word from the PDP-11. The threshold voltage is manually variable at 
present but a simple modification can place this under computer control 
also. 

The horizontal and vertical blanking pulses reset all counters in the 
thresholder, and phase lock its master clock. The DATAIN flag shown in 
Fig. 5.2„1 i s set to l whenever a data word: is ready for transfer to the , 
frame.* generator. DATAIN i s rese t to 0 by the frame generator upon receipt 
of the transferred data word. Only 15 TTL packages are required to imple- 
ment this thresholder. 



Figure 5.2.1 Diagram of Vidicon Interface 
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•5. 3 DIGITIZED FRAME GENERATOR 

The digitized frame generator shown in Fig. 5. -3.1 forms the final digi- 
tal picture cells by combining data from successive horizontal scan lines 
as shown. Although 1 TV frame is made up of 2 fields (2:1 interlace) 
only 1 field is needed in this process since the field, resolution is 
only 250 lines . 

The control flow of the frame generator is shown in Fig. 5.3.2 . Be- 
fore entering the main control loop all resolution information is loaded 
from the PDP-11, and all. counters are initialized. Within the control 

\ 

loop the generator processes 

1. Data words from the video thresholder, forming the picture elements. 

•2. Control words from the PDP-11, and 

3. All counters and flags during the horizontal blanking pulse. 

The frame generator is constructed exclusively from Register Trans- 
fer Modules (RTMs) . These modules allow hardware to be designed at the 
register transfer level; that is, the level at which data words (16 bits) 
are transferred between storage locations and operated upon. Table 1 
lists the RTM data modules used to construct the frame generator. It also 
serves to define all registers and flags. Figure 5.3.3 shows the RTU imple- 
mentation of the main control flow in Figure 5.3.2. Although Figure 5.3.3 
resembles a software flowchart it has, in fact, been directly implemented 
in hardware via backplane wiring between the various RTM modules. 

To provide some tutorial background on RTMs, the three types of 
modules used in RTM implementations are listed below with brief functional 
descriptions. 


- — •' — ' — 

bights and Switches 

■ -■ — 

displays bus transfers 

Bus Control 

reg : BSR 

controls bus use. 

General Purpose 
Arithmetic Unit 

reg: A,B 

performs arithmetic 
operations 

16 word scratch- 
pad nicmdry 

V 

• V *• 

f 

. 

reg: 

IN 

picture element pointer 

""LINTOT 

horizontal line counter 

LIME 

picture element line counter 

L1MEMAX 

number of lines in picture 
element. 

IK scratch pad 
memory 

reg: 

MA 

memory address 

"MU 

memory buffer 

4 word read only 
memory 

reg: 

MASK 

masks PDP-11 control word 

2<0 10 

- 

number of horizontal scan 
lines in 1 video field. 

Boolean Flag 
Modules 

i 

i 

• • 

■: • 

• 

• 

flag: 

infl ag 

PDP-11 control word ready 

outf lag 

output word for PDP-11 ready 

firstline 

horizontal line is 1st line 
of picture element 

start 

s tarts/stops digitized 
frame generator. 

skipline 

set to skip 1st line of 
video field. 

lastline 

horizontal line is last 
line of picture element. 

first frame 

used to skip every other 
video field. 

linedone 

disables counter updates 
during horizontal blank time. 

Interface 1 

reg: 

IN REG 1 

input from PDP-11 

0UTREG1 

output to PDP-11 

Interface 2 

reg: 

IN REG 2 

input from video thresholder 

OUT REG 2 

output to video thresholder 

. Video Thresholder 

flag: 

— 

SAMPLE 

data word ready 

NEW FRAME . 

positive vertical blanking pulse 

NEWLINE 

positive horizontal blanking 
pulse. 


Table 1 . RTH Modules for Digitized Frane Generator. 
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1. Control modules - these evoke data operations in their proper sequence, 
and select paths through diverging control sequences on the basis of 

l . 

boolean flags. 

2. Data modules - these perform typical data operations (addition, shift- 
ing, boolean logic) on data words. 

3. Memory modules - these are used for data word storage. 

As a further illustration of the RTM concept, the details of the 
Kmacro (SAMPLE) operation are shown in Fig. 5 .3 .4 .This portion of the 
control flow is executed when a data word from the video thresholder 

\ 

buffer is to be processed. The work is OR’ed with data from previous 
horizontal scan lines to form the picture .cells . 
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5,4 PDP-11 PROCESSING 

The PDP-rll minicomputer must be able to accept the picture cells 

V * 

from the digitized frame generator, search them for Is . (points of light) 
and record the location of these light intensity peaks. This data can 
be displayed directly on a teletype or graphics display terminal or in 
some statistical fashion. Eventually the PDP-11 will be used as a real- 
time controller for the optical system, with direct control of the input 
image, spatial filter, and digitized frame generator. 
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CHAPTER 6 

EXPERIMENTAL RESULTS 

6.0 Introduction 

As has been pointed out earlier in this report, no real-time data from 
phased array antennas was made available to us for a variety of reasons. It 
was thus reluctantly decided to simulate the phased array data. To this end 
a special purpose function generator was constructed to produce simulated input 
plane patterns. This section describes the simulator and provides documenta- 
tion on the real-time processing of radar data using the Pockels tube input 
modulator (Chapter 1), the input formats described (Chapters 3 and 4) and the 
computer interface (Chapter 5). 

6 . 1 Linear Phased Array Simulator 

We recall from paragraph 3.3 that the returned signal at the nth element 
of a linear array, after heterodyning, has the form 


v in (t) = P T (t) cos + f d > e - 2™“] 

; P 


where the time origin t = o is taken at T = 2 R/C, the round trip to the 

target, f^ is the intermediate frequency of heterodyning, f^ is the Doppler 

frequency shift, and a = - 4 - sin 8 is the differential delay per element in 

’ c 

cycles of f c (carrier frequency) and d is the uniform spacing of the array 
elements, X is the RF carrier wavelength, 0 is the target angle from ore- 
sight and Tp is the pulse duration. 

The radar simulator produces these waveforms as a train of square pulses 
at a frequency (f^+ with a total duration T^ and a linear phase shift per 
line. Tp is adjustable in the range 0 to 50 ^sec, while the frequency of the 
pulse train can vary from 20 KHz to 10 MH 3 . The radar simulator is phase-locked 


r-~ 


SYNC. GENERATOR 



Fig. 6.1.1 Black Diagram of Linear Phased Array Simulator 
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to a timing signal to allow any variable delay per line At to be introduced 
while maintaining the proper phase relationships between signals on differ- 
ent lines. (Recall that the signal written on each line of the electro-optic 
target crystal corresponds to one of the antenna array elements.) 

Linear arrays of up to 256 receiving elements can be simulated (with 
extension to 500 elements possible) ; the number of elements and their spacing 
can be varied by control knobs on the front panel of the simulator. The 
position of the returned pulse signal in the horizantal time axis, as well as 
the pulse duration and the relative delay between consecutive array elements, 
can all be independently varied by control knobs. 

The functional diagram of this phased array radar simulator is shown in 
Figure 6.1.1 

6.2 Basic Optical System 

Figure 6.2.1 shows the basic optical system used. To conserve table 
top space, a stand for the laser is suspended from the bottom of the table 
top. The laser, a focusing lens and aperture are mounted on this stand. 

The remaining optical system is on the table top. A periscope prism arrangement 
brings the beam to the table top. A fairly plane, uniform intensity beam over the 
50 x 50 mm dimensions of the target crystal is needed for optical data pro- 
cessing. This 70 mm diameter (50/2) and the laser beams size to the 1/e points 
determines the collimating optics. A 20X AR coated microscope objective focuses 

the beam to a diffraction limited spot. A spatial filter pinhole was found to 
be essential if a uniform beam was to be obtained. A 10pm diameter pinhole works 
well to eliminate noise in the beam without loosing laser power. The incident 
beam diameter is 2 mm, the 20X objective and a 381 mm collimating lens provide 
a 75 mm collimating lens provide a 75 mm diameter beam which illuminates the 
crystal. An additional input polarizer is used to insure linearly polarized 
incident light. The output polarizer is crossed to provide no light out with 
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no voltage on the crystal. 

The spatially modulated collimated laser beam is now Fourier transformed. 
The focal length of the transform lens is determined by matching the optical 
spatial bandwidth to that of the vidicon used for pick-up and detection. The 
spatial frequency n is related to the focal length f by 


n = x/Xf 


where x is the distance in the tramsform plane associated with the vidicon 
scan, X is the wavelength of the laser light (633 nm) and f is the focal 
length of the Fourier transform lens. The electro-optic crystal resolution 
is 20 lines/mm and the vidicon scan dimensions are 0.8 cm x 0.6 cm, from which 
the focal lengths needed are 


f, . < 635 mm 

honz - 


f ^ < 475 mm 
vert — 


These limits insure an optical spatial bandwidth greater than or equal to 
the highest spatial frequency on the Crystal. With the 495 mm focal length 
transform lens used and a 2.5 cm diameter vidicon, a resolution of 


^oriz = ^5*6 lines/mm 

n ^ = 19.2 lines /mm 
vert 

can be resolved in the spatial frequency plane. These values are comparable 
with the vidicon and crystal specifications. 

The Fourier Transform plane can be reviewed directly on the vidicon or 
first magnified and then viewed on the vidicon (in the case where the maximum 
spatial frequencies contained on the input pattern are lower than the maximum 
spatial bandwidth allowed by the crystal). In the second case, a scanning 
microscope is used which can also drive an X-Y platter for beam recording. With 


. 1 *' 
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the use of the magnifying system, enlarged output plane patterns can be pro- 
jected on a screen and pictures taken while the vidicon drives the computer 
interface and display. 

6.3 Results 

On line simulated radar returns for a linear array have been written on 
the input KD*P crystal and the corresponding output plane formats analyzed. 

The radar signal simulator described in paragraph 6.1 has been used through- 
out the experiments. The optical configuration used was shown in Fig. 3.3.1. 

The input collimated plane waveform is modulated by the electrical signal from 
the radar simulator which is written on the electro-optic crystal. With this 
crystal between crossed polarizers an on-line two-dimensional amplitude 
modulated light beam results. The optical system forms the Fourier transform 
of the image on the crystal and displays a magnified version of this transform 
on the TV monitor. The RTM computer interface then digitizes this output plane 
pattern on line with a variable cell size and threshold level. 

All results correspond to a linear phased array of 100 elements spaced 
at a distance equal to one half the RF carrier wavelength. A pulsed RF signal 
of 20 ysec duration. is used in all cases. The intermediate frequencies used 
are 1 MHz and 500 KHz. 

Figures 6.3.1 to 6.3.8 show the results of this initial on-line processing 
at various points in the system and for several input targets. All of the pro- 
cessing and data shown are on-line. 

Figure 6.3.1 shows the output plane pattern for a boresight target and a 
radar intermediate frequency of 1 MHz. Note that the first order diffraction 
peaks lie on the vertical axis, (u axis) 

Figure 6.3.2 shows the output plane pattern for a target 7° off the bore- 
sight with an f ^ = 1 MHz. Note that the location of the first order diffraction 


peaks are now displaced off axis along the v direction proportional to sin(7°) . 

Figure 6.3.3 is the RTM digitized version of Fig. 6.3.2 with suppressed 
D.C. term. This prevents saturation of the vidicon pick-up tube used. 

Figure 6.3.4 is the output plane pattern for a target 20° off boresight 
with an f = 1 MHz. Note that the shift in the first-order diffraction peaks 
is larger than before, (proportional to sin(20°)). 

Figure 6.3.5 is the RTM digitized version of Fig. 6.3.4 with the D.C. 
term again suppressed. 

Figure 6.3.6 is the output plane pattern for a target at 10° off bore- 
sight and a different radar IF frequency f^ = 500 KHz. 

Figure 6.3.7 is the RTM digitized version of Fig. 6.3.6 with the D.C. 
term only partially suppressed. 

Finally, Fig. 6.3.8 corresponds to the same input pattern than Fig. 6.3.2 
but using the phase-modulation capability of the electrooptic crystal using 
only an input polarizer. 

These initial results prove the on-line beam forming capability of our 
KDP light modulator in radar processing. The validity of the input format 
described in the prior reports and the output plane pattern predicted by this 
theory were also confirmed by these tests. The on-line operation of the RTM 
computer interface in an on-line radar data situation has also been demonstrated. 
The output plane information is also available as coordinate pair locations of 
all principal maximums. From such data, the computer can automatically extract 
target angle information (from the v displacement of the spots) ; this data can 
then be used to further control the radar system or other supporting radar 
systems, for data storage or for display purposes as demonstrated. 
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CHAPTER 7 
FUTURE RESEARCH 

7.0 Introduction 

In our research on radar processing using an electro-optical 
processor we have also considered the possibility of processing 
linear FM stepped, pulsed Doppler, and CW Doppler returns. The follow- 
ing three sections briefly outline the procedure by which the electro- 
optical processing techniques described in the report could be applied 
to this alternate radar systems. 

The remaining sections of this Chapter outline more advanced 
processing techniques for radar data using standard optical processing 
techniques such as matched filtering, spatial filtering and correlation. 
One of the key elements in this advanced applications of optical process- 
ing to radar data involves a survey and experimental feasibility efforts 
in the area of on-line spatial filter materials. Section 7.5 discusses 
one seemingly viable material for such an adaptive on-line spatial filter. 
This Chapter and the aforementioned discussions will purposely be kept 
very brief as they are quite numerous and constitute our proposed research 
program in optical processing of radar data for the coming year. 

7.1 OPTICAL PROCESSING OF PULSED-DOPPLER DATA 

In pulsed Doppler systems a set of single frequency pulses is 
transmitted at regular time intervals. For a constant velocity target 
in a noise-free environment a linear phase slope will again exist across 
the returned echoes as it was the case in the phased array processing 
described in Chapters 2, 3, and 4. With a pulsed-Doppler system the 
target’s velocity can be determined from this linear phase slope. 
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7.1.1 PULSED DOPPLER PROCESSING IN THE PRESENCE OF NOISE 

In the presence of noise the phase slope across the returned 
echoes can be optically determined by correlating the input data with 
the phase slopes corresponding to all possible velocities within a 
given range and resolution bin. Fig. 7.1.1 shows the expected correlation 
peak provides information on the target velocity. 

7.1.2 OPTICAL PROCES SOR 

Using our Pockels tube modulator (Section 6.2) and the conventional 
optical processing configuration in Fig. 3.3.1 the returns for consecutive 

v 

pulses will again be written on separate raster lines on the input crystal 
after being heterodyned and range gated . The two dimensional optical 
Fourier Transform of this input pattern as before has a first order 
maximum at back plane coordinates from which information on the target’s 
velocity can be obtained. 

Since N resolution cells of target velocity exist between the grating 
lobes of the transform it is possible to simultaneously on-line in- 
parallel correlate against N different velocities in real time. Multi- 
target velocity resolution within range bins also appears possible. 

The computer interface plays a major role in such a system since 
one might envision a totally automated processor in which one first obtains 
coarse resolution by one cycle through the system. This coarse resolution 
information could then be used to determine the fine range spatial filter 
required to be inserted in the spatial filter plane. These filters may 
be computer generated ones which are stored in the computer memory. A 
more practical approach within the scope of current technology would 
involve storing these filters optically as multiple spatial filters 
(holograms). This will be expanded in Section 7. 
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1.2 QP-TTCAL PROCESSING OF LINEAR-STEP FM RADAR DATA 

A second type of radar signal which we considered processing was 
Linear stepped FM data. In this radar mode M consecutive RF pulses 
are transmitted at successively higher frequencies. A. wide frequency 
spectrum is thus covered. Again a linear phase slope will exist across 
the returned echoes in this system for stationary targets. This slope 
is now a function of the target’s position within a range resolution cell. 
7.2.1 OPTICAL PROCESSOR 

Using a similar optical configuration as before and a similar input 
format the returned echoes can again be correlated with corresponding 
phase slopes. The Fourier transform of the input pattern described 
has a first order maximum from which one can now obtain fine range 
resolution. With N resolution cells between continuous grating lobes, 
one can again perform a parallel optical correlation against filters of 
N ranges and thus obtain fine range resolution in any. of N range resolution 
bins. 

The correlation filters, the optical spatial filters and the role 
of the computer interface are various other considerations and factors 
in this processing that must be investigated. 

7.3 OPTICAL PROCESSING OF CW DOPPLER DATA 

Optical processing of CW doppler data directly parallels the 
pulsed doppler scheme. With continuous data rather than pulsed wave- 
forms, the input processor is even simpler. 

7.4 REAL-TIME DATA 

Until actual phased array radar return data is available, efforts 
should be directed at the aforementioned schemes since real-time data 
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in these forms exists. The questions that remain to be answered are still 
quite numerous. 

(1) The degree of preprocessing and post processing required varies 
With the radar type and system and must be separately evaluated for 
each case. 

(2) The optimum optical system with special lenses etc, must be 
determined. This again varies with the radar system chosen. 

(3) The spatial filter resolution necessary must also be determined. 
The processing times required will also affect the spatial filter material 
chosen. 

(4) The resolution, speed, number of gray levels required, etc. 
in the output plane and the data desired from the computer interface 
should be determined for each radar processing made envisioned. Extension 
and expansion of the interface will then be possible. 

7.5 SPATIAL FILTER MATERIALS 

» 

Many of the radar processing applications aluded to above involve 
spatial filtering, matched filtering,- and correlations. It is possible 
to perform these powerful operations in parallel in an optical 
processor. Among the more widely researched materials for use in optical 
storage of information are: 

1. Film 

2. Photothermoplastics 

3. Liquid Crystals 

4. Photodichroic crystals 

Film has the disadvantage of not being adaptive and of not being 
on-line. Photothermoplastics generally have a lifetime problem 
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that limits them to 10,000-100,000 cycles, liquid crystals are generally 
adequate for display purposes but their slow response times are not com- 
patible with the on-line optical processing envisioned. Photothermoplastics 
generally require larger optical power levels but indications are that 
this problem can be reduced, fatigue problems overcome, etc. We are 
currently investigating them in an alternate application but feel they 
deserve attention in this specific contest of radar processing. 

7 . 6 COMPUTER INTERACTIONS 

In perhaps the most all encompassing project, the returned signals , 
from a phased array, linear FM stepped, and pulsed doppler radar systems can 
be processed on one system with the digitized output plane, the computer 
interface, and the PDP-11 being used to scan the results on one radar 
display and store it and use the result to control a linear FM stepped 
radar to obtain range information, then a pulsed doppler system to obtain 
velocity data. The interaction of such a number of systems may be easier 
than at first appears. Use of the computer to store a fixed program cor- 
responding to the returns from such a set of 2 or 3 radar systems seems 
reasonable. Obtaining actual radar data from 2 or 3 systems for the same 
target seems an unreasonable request thus simulated returns would be used. 

The scope of such an effort and many details are clearly rather involved — 
but perhaps a reasonable attempt at such a program could be made with 
initial results from the PDP-11 feeding our modulator expected in one year. 

In a more limited and controllable situation one could use the output 
from the PDP-11 to control the fine and coarse range controls noted in 7.1.1. 

By using the digital processor to control the updatings of the spatial filters, 
the input data, and the stored and transmitted output data, one can 


- 109 - 


demonstrate the computer feedback control within a well defined and 
controllable and yet useful context. 

Such a system provides a digital optical radar processor much less 
complex and cheaper than existing alternatives such as a dedicated lai*ge 
scale computer, should greatly advance the state of the art in coherent 
optical data processing with on-line spatial light modulators. 


SUMMARY 


An on-line spatial light modulator developed at CMU suitable as the input 
transducer for a real-time optical data processing system has been described 
(Chapter 1). 1000 x 1000 point resolution and on-line operation at TV frame 

rates is possible. The use of such a device in the analysis and processing of 
radar data in real time has been considered' and initial results reported. 

The fundamentals of signal processing for linear (Chapter 3) and planar 
(Chapter 4) phased array radar data as well as Pulsed Doppler and Linear-Step 
FM system (Chapter 7) using this input modulator have been developed. Tima 
and space multiplexing input formats, full and thinned aperture processing, 
and amplitude and phase modulation schemes using the system are developed, 
analyzed and compared. 

An interface from the optical processor to a control-digital computer 
has been designed, constructed, and tested. The ability to digitize the entire 
Fourier transform plane or the correlation plane of an optical processor in 
real time and input the resultant data to a digital computer for storage or 
for upe in controlling the entire optical processor has been demonstrated. 

The cell resolution and size and the gray scale levels in the backplane are 
variable under program control. 

Real-time processing of linear phased array under data has been performed 
for the first time by other than an acousto-optic modulator. The Pockels 
tube modulator offers advantages in fabrication, angular coverage and resolu- 
tion, in the larger array size it is now possible to handle and in-depth of 
modulation etc. 

The input transducer, optical system, and computer interface have all been 
operated in real-time with real-time radar data with the input data returns 
recorded on the input crystal, processed by the optical system and the output 
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plane pattern digitized, thresholded and outputted to a display and storage 
in the computer memory all in real-time synchronized to TV frame rates. 

No other input modulator offers the real-time operating characteristics 
of the Pockels tube. Its ability to process radar data in real-time has been 
demonstrated.. The correlation with theory and experimental results have been 
most encouraging.. Recommended directions for future research in the area of 
optical processing of radar data in real-time are sketched in Chapter 7 and the 
proposal. 
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